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We theoretically investigate collision of optical beams travelling in opposite 
directions through a centrosymmetric photorefractive crystal biased by a 
spatially non-uniform voltage. We analytically predict the fusion of coun- 
terpropagating solitons in conditions in which the applied voltage is rapidly 
modulated along the propagation axis, so that self-bending is suppressed by 
the "restoring symmetry" mechanism. Moreover, when the applied voltage is 
slowly modulated, we predict that the modified self-bending allows conditions 
in which the two beams fuse together, forming a curved light-channel splice. 
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Photorefraction has been a principal setting for the experimental investigation of soliton 
collisions [1]. One setback in previously reported schemes is the fact that the single solitons 
intrinsically self-bend along their propagation [2]. Although this has a minimal effect on 
spatial soliton interaction where both beams propagate in the same direction, it strongly 
affects the phenomenology when the beams counterpropagate [3,4]. For example, in soliton 
head-on collisions, the main physical effects emerge when there is an extended beam spatial 
overlap [5-7], and this overlap is strongly reduced by self-bending that curves both beams in 
the same lateral direction. A possible solution to enhance spatial beam overlap is to use spe- 
cific curved trajectories obtained by launching tilted beams [4]. Both from the fundamental 
and applicative perspective, it is natural to formulate an alternative strategy allowing the 
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increased overlap without a specific tilted alignment. For example, a self-adaptive fusion of 
two counterpropagating beams could amount to an alignment-robust splicing technology for 
fiber connections. 

Recent studies have shown how self-bending can in fact be suppressed by the "restoring 
symmetry" mechanism implemented through the use of alternating bias voltage profiles along 
the propagation direction [8,9]. In this Letter we theoretically consider the head-on collision 
of two counterpropagating optical beams travelling through a centrosymmetric photorefrac- 
tive crystal in the presence of a system of electrodes delivering a voltage profile periodically 
modulated along the propagation axis. When the period of the applied voltage is much 
smaller than the optical diffraction length, the consequent suppression of the self-bending of 
each single beam allows us to investigate the interaction of two micron-sized counterprop- 
agating beams whose complete spatial overlap, along a straight line, triggers their mutual 
longitudinally nonlocal interaction [4,5]. In the situation where the two beams are mutually 
incoherent, we show that soliton fusion occurs by analytically proving the existence of a two- 
parameter family of fully overlapping counterpropagating solitons. In conditions where the 
applied voltage is slowly modulated, we numerically identify the conditions for exploiting the 
longitudinal wiggling beam profile [9] to achieve a robust fusion of two counterpropagating 
beams impinging on the opposite crystal facets. The results indicate a self-adaptive merging 
of two channels, an effect that has a potential photonic application in minimizing optical 
power loss associated to diffraction and misalignment in optical fiber splicing. 

Consider a photorefractive crystal whose facets x = —L x and x = L x have a set of 
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electrodes that deliver the periodical potential profiles — V cos(k v x) and V cos(k v x), re- 
spectively, as shown in Fig.(l). In Ref. [9], it has been proved that, if I(x,z) is the optical 
intensity of the light travelling along the z-axis through the crystal, the photorefractively 
induced refractive index change is 
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where x — K B T/q, a = — (l/2)nlgel(e r — l) 2 (T is the crystal temperature, q is the electron 
charge, e r is the relative dielectric constant at the given T, uq is the uniform refractive 
index background, g is the significant quadratic electro-optic coefficient), is the intensity 
of a reference background uniform illumination, and ip = VoIb/[L x cosh(K v L x )]. In the TE 
configuration, the complex amplitude of the monochromatic (at frequency to) optical electric 
field E(x,z) satisfies the Helmholtz equation (d xx + d zz )E + kl(n + 5n) 2 E = where 
ko = uj/c and 5n is given by Eq.([T]). In order to describe head-on collision of two beams 
counterpropagating along the z-axis we set E(x, z) = exp(ikz)A + (x, z) + exp(— ikz)A_(x, z) 
where k = kono and A + and A_ are the slowly-varying amplitudes of the forward and 
backward propagating beams, respectively. Inserting this expression for E into the Helmholtz 
equation, in the paraxial approximation and noting that for k v « k (i.e. the period of 
voltage modulation is much greater than the optical wavelength) light cannot be Bragg- 
matched with the periodic refractive index profile, A± satisfy the coupled parabolic equations 
[±id z + (l/2k)d xx ]A± = —{k/no)8nA±. We here focus our attention on the relevant case of 
two mutually incoherent counterpropagating beams for which the total optical intensity is 



given by I = \A + \ 2 + | ^4 | 2 . Since Eq.flT]) is not an even function of x if the intensity / is 

even, the two beams A± do not propagate on a straight line: they experience the effect of 
self-bending in the same lateral direction, so that the effect of their interaction is generally 
limited by the smallness of the overlapping region. In order to maximize the overlap of the 
two beams we take k v 3> 2n/L d (the situation corresponding to the geometry illustrated in 
Fig. (la)) where Ld is the longitudinal scale characterizing the propagation of A± (Ld typically 
coincides with the optical diffraction length). In these conditions the optical beams are not 
able to follow the rapid voltage oscillation and the averaged fields do not experience self- 
bending, and do not wiggle (the mechanism of "restoring symmetry" discussed in Ref. [9]). 
In this regime (i.e. k v ^> 2ix / L d ) it is possible to set A± = ^fhV± + 5A±, where V± are those 
parts of the fields having a longitudinal scale of variation Ld, and 5A± are longitudinally 
rapidly varying, on a scale 2tt/k v , and they are uniformly in the condition \5A±\ <C \/7b|V±|. 
This self-consistent decomposition of the fields into a slowly varying mean-field component 
and a rapidly oscillating and small correction allows us to derive a set of equations for V± 
(following a procedure very close to that reported in Ref. [9]) that are 
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where we have also introduced dimensionless variables according to £ = k\i/j/If,\y2\a\/noX, 
( = k(\a\/no)(^/Ib) 2 z and 7 = 2/c 2 x 2 |a|/ri - Note that, as expected, the two beams V± 



are driven by the very same nonlinear waveguide (as a result of the two beams mutual 
incoherence) which, if \V + \ 2 and \V-\ 2 are transversally even functions of £, is transversally 
even as well so that no self-bending occurs. Equations admit of the solution 



for any values of the real parameters $ and a if the function v(r) satisfies the equation 



medium in the presence of the "restoring symmetry" mechanism, as discussed in Ref. [9], so 
that the fields in Eqs.Q constitute a two-parameter family of counterpropagating solitons. 
It is worth stressing that the two solitons of each pair do not suffer self-bending, are fully 
overlapping and therefore Eqs.(j3]) describe fusion of solitons counterpropagating along a 
straight line. The parameter $ sets the mutual power content of the two solitons in such a 
way that (V^l 2 + \ V-\ 2 = \v\ 2 , whereas the parameter a (subjected to the restriction a C 1 
required by the paraxial approximation) allows the soliton pair to be slightly tilted with 
respect to the z-axis. 

In order to check the above analytical results and to extend our investigation to the off 
axis interaction configuration (see Fig. (lb)), we have integrated the full time-dependent pho- 
torefractive nonlinear optical model [10]. In our numerical approach, at each instant of time, 
we evaluate the electric field distribution induced by the boundary applied voltage and the 
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Note that Eq.fllJ) coincides with the equation describing solitons propagating through the 
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photoinduced charge solving the (x, z) electro-static Poisson equation, and the correspond- 
ing optical field distribution determined by the electro-optic response through the parabolic 
equation [11]. We have chosen a crystal bulk (layer) of potassium lithium tantalate niobate 
(KLTN) (no = 2.4) of thickness 2L X = 2 x 50 fim and length L z = 1000 um. In order to 
investigate fusion of coaxial counterpropagating beams in the fast modulated regime (with 
electrode modulation period 2tt/k v = 200 um), we have chosen 7 = 0.2 (an experimentally 
available situation as reported in Ref. [12]) and we have launched two identical counterprop- 
agating Gaussian beams A+(x, 0) = v/l&/2/(£) at z = and L) = ^Ib/2f(£) at z = L 
where /(£) is a real Gaussian profile centered at £ = (£ is the same dimensionless spatial 
variable as in Eqs.(f2"l)). We have performed various numerical simulations varying both Gaus- 
sian width and amplitude together with the applied voltage thus determining their values 
which maximize the overlap between the forward and backward propagating field profile at 
z = 0, so to observe the formation of a stable and straight optical channel (see Fig.(la)). 
In Fig. (2a) we have plotted the FWHM, a, of / 2 (£) as a function of / = /(0) (plotted 
as stars) corresponding to coaxial beam fusion and, for comparison purposes, we have also 
reported the theoretical soliton existence curve (solid line) derived by Eq.flU) (see Ref. [9]). 
The good qualitative agreement between the two different situations indicates that the fusion 
mechanism is robust and feasible. 

If the applied voltage is slowly modulated, a form of modified self bending occurs since the 
optical beams are able to adiabatically follow the electrode modulation [8,9]. This property 
can be profitably exploited to design a configuration where two beams, impinging onto the 

7 



crystal facets z = and z = L along two parallel propagation directions, are made to form a 
single and curved light channel within the crystal bulk. In the circumstance of the geometry 
depicted in Fig. (lb), for example, the applied voltage is reversed one time along the z— axis 
so that beams fusion is possible since the optical beams bend toward negative and positive x 
direction in the z < 500/im and z > 500/U.m crystal regions, respectively. We have investigated 
off axis beam fusion in the configuration as in Fig. (lb) by means of the above discussed 
numerical scheme by launching two identical but shifted Gaussian profiles (i.e. by setting 



A+(x,0) = A exp[{x - d/2) 2 )/{2s 2 )} at z = and A-(x,L) = A exp[(x + d/2) 2 )/(2s 2 )] 



at z = L, where A = 2^/4 and s = 3 fj,m) for various different applied voltages V and 

mutual beam displacement d thus determining their values which maximize the overlap 
between the forward and backward propagating field profile at z — 0. The result of these 
calculations are reported in Fig. (2b) from which we note that fusion can be attained even 
for counterpropagating beams whose mutual distance d is much greater than their common 
width. This result suggest a feasible way for obtaining self-adaptive optical fiber splicing. 
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List of Figure Captions 

• Figure 1: Geometry of the collision between counterpropagating optical beams (re- 
ported as shaded regions around x = 0) through a non-uniformly biased photorefrac- 
tive crystal layer (black and gray stripes are here electrodes at opposite potentials), 
(la) Fusion of two coaxial counterpropagating solitons in the fast modulated regime, 
(lb) Merging of two off axis counterpropagating beams into a single optical channel 
due to modified self-bending in the slowly modulated regime. 

• Figure 2: (2a) Intensity full width at half maximum FWHM a as a function of the peak 
amplitude fo of the Gaussian input counterpropagating beam profiles allowing beam 
fusion (stars) and corresponding theoretical existence curve (solid line) associated with 
counterpropagating soliton fusion (evaluated from Eq.(j4])). (2b) Voltage V as a function 
of the displacement d between the two beams required to form an optimal fused splice 
along a curved trajectory. 
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Fig. 1. Geometry of the collision between counterpropagating optical beams (reported as 
shaded regions around x = 0) through a non-uniformly biased photorefractive crystal layer 
(black and gray stripes are here electrodes at opposite potentials), (la) Fusion of two coax- 
ial counterpropagating solitons in the fast modulated regime, (lb) Merging of two off axis 
counterpropagating beams into a single optical channel due to modified self-bending in the 
slowly modulated regime. 




Fig. 2. (2a) Intensity full width at half maximum FWHM a as a function of the peak ampli- 
tude f of the Gaussian input counterpropagating beam profiles allowing beam fusion (stars) 

and corresponding theoretical existence curve (solid line) associated with counterpropagat- 
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ing soliton fusion (evaluated from Eq.(j4])). (2b) Voltage V as a function of the displacement 
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